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Abstract: A 1 /L-expansion for percolation problems is proposed, where 
L is the lattice finite length. The square lattice with 27 different sizes L = 
18, 22, . . . 1594 is considered. Certain spanning probabilities were determined 
by Monte Carlo simulations, as continuous functions of the site occupation 
probability p. We estimate the critical threshold p c by applying the quoted 
expansion to these data. Also, the universal spanning probability at p c for 
an annulus with aspect ratio r = 1/2 is estimated as C = 0.876657(45). 

1 Introduction 

In reference [T] a square lattice is viewed as a torus, thus without frontiers, 
where sites are randomly occupied with probability p. Wrapping percolation 
probabilities can be defined within this geometry, counting configurations 
which wrap along the horizontal and/or/xor vertical directions. The proba- 
bility Rh, for instance, counts all configurations wrapping along the horizon- 
tal direction, no matter which occurs vertically. As a function of p, Rh cor- 
responds to a plot like figure 1 (to be quoted later) for a finite lattice. In the 
thermodynamic limit L — > oo, this plot approaches a step function, Rh = 
below, and Rh = 1 above the critical threshold p c . The point Pl{ t ) exem- 
plified on the plot serves as an estimator for p c : by measuring a sequence of 
such values, for larger and larger sizes L\, L 2 . . . Ljv, one can extrapolate this 
sequence for L — ► oo. Reference [T] presents the figure p c = 0.59274621(13), 
the most accurate available today, obtained from more than 7 x 10 9 Monte 
Carlo samples. Starting from an empty lattice, filled site by site at random, 
each sample provides a single number to the statistics, namely the precise 
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number n of occupied sites for which the proper wrap (horizontal for Rh) ap- 
pears for the first time. This computational strategy of filling up the lattice 
site by site and storing data onto bell-shaped n-histograms is equivalent to 
that of 2j, and similar to early works [3 . However, the multi-step strategy 
introduced in pP involves many other components (see |3]). The authors of 
[T] assume a Pl( t ) — Pc ~ L _ ( 2+1 /^ dependence, where v = 4/3 is the corre- 
lation length critical exponent. As 2 + 1/v = 2.75 is a sufficiently high value, 
data up to only L = 128 were needed. Also, the authors adopt the value 
r* = 0.521058290, corresponding to that particular wrapping geometry and 
exactly known [5] as the limiting value for Rh at p c , called Pinson's number in 
0. Outside this value or within other geometries, the not-so-high exponent 
1 + 1/u or even the smaller 1/V were observed [7] . 



Figure 1: Spanning probability function for L = 18 (solid line), 38 and 74 
(dotted lines). By fixing some value r, one can find a sequence of values 
Pl{t) for increasing lattice sizes, approaching the critical threshold p c . 

Here, we propose the mathematical form 
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for estimators pi obtained from quantities like Rh, where the cutoff M is 
conveniently chosen according to the numerical accuracy available. 
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Table 1: Values for p L obtained for fixed r = 0.9, as an example. A Chi- 
square fitting of equation (1) gives p c = 0.59274675(88), A = -0.44204(11), 
Ai = 3.275(14), etc. 

We apply this formula to the L x L square lattice, with a set of different 
increasing lengths such that the numbers of sites grow by a factor of \/2. 
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Table I shows an example for r = 0.9, for which the traditional Chi-square 
fitting [H] gives p c = 0.59274675(88), in agreement with P although within 
a larger error bar. We adopted M = 4 in equation (1), compatible with 
our smallest lattice size L = 18, since 18~ 4 ' 75 = 1 x 10~ 6 still falls inside 
our numerical accuracy, whereas the next term 18~ 5 ' 75 = 6 x 10~ 8 would be 
outside. The quality of this fit can be appreciated by the so-called goodness- 
of-fit Q |E], a quantity between and 1. The fit is considered believable jS] 
for values of Q > 0.1. In our example, we get Q = 0.857 for table I. All 
our data to be discussed hereafter, for many other values of r between 0.5 
and 0.99, present the same degree of accuracy, giving credit to our proposal, 
equation (1). In spite of these accurate results, one cannot rule out some 
possible higher terms deviating from (1), for example that proposed in jUj. 

2 Measured quantity 

First, let's explain the spanning probability we adopted within the torus, 
instead of the wrapping probability pQ. We consider two parallel lines distant 
L/2 from each other, on the L x L square lattice. For each sample — again 
obtained by filling up the initially empty lattice, site by site at random 
- we count the precise number n of occupied sites for which these lines 
become connected for the first time, no matter which occurs around the other 
direction. This approach has a big advantage over the wrapping probability 
around the whole torus: From the same sample we can count n just L times 
instead of only once! The parallel lines can be numbered + L/2) for 
% = 1,2, ... L/2 along the horizontal direction, with the same procedure 
repeated vertically. Thus, the statistics is multiplied by a factor of L. In 
table I, for instance, the sampling counting 10 9 for L — 18 corresponds to an 
n-histogram with 18 x 10 9 accumulated units (the total area below the bell- 
shaped curve). In the same table, for L = 1594, the much smaller sampling 
counting 4 x 10 6 corresponds indeed to almost the same statistics, i.e. 6 x 10 9 
accumulated units below the curve. This trick allows us to test a wide range 
of lattice sizes, and verify the validity of our proposal (1). The further 
computational time one needs in order to implement this trick is negligible: 
we simply keep in memory the top and bottom (left and right) extreme lines 
for each already formed cluster of neighbouring occupied sites. Thus, for 
each new included site, only the last updated cluster must be verified. 

The spanning probability function is obtained by superimposing a lot 
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of step functions, one for each counted n, and dividing the result by the 
total number (Lx sampling counts). An average is then performed, weighted 
by C(L 2 ,n)p n (l — p) L ~ n where C is a combinatorial factor, yielding the 
p-continuous curves shown in figure 1. 

The 3-digits error bars shown in table I were obtained by dividing the 
whole set of data into, say, S — 10 sub-sets, independently calculating Pl(t) 
for each sub-set. The error bars are the standard deviation of this distribution 
divided by y/S. This last division is based on the supposition that the whole 
data is normally distributed. In order to verify the validity of this approach, 
we repeated the same procedure with S = 20. Indeed, the error bars are 
approximately the same, independent of S. For safety, we adopted always 
the largest between both error bars so obtained. For intermediate lattice 
sizes (from L = 74 up to 402), we used S = 10 and 5, instead of 20 and 10. 

We also simulated larger L x L lattices for r = 0.5, up to L = 24000, 
with free instead of periodic boundary conditions, within a poorer statistics. 
The results are also compatible with equation (1). 

3 Data Analysis 

By fitting data with equation (1), we have a set of M + 2 parameters to 
be determined, namely p c , A (t), Ai(t), . . . Am(t). The error bars for these 
quantities come from the Chi-square fitting procedure [S] , as a consequence of 
the primary error bars directly measured for the crude data (that of table I, 
for instance). Thus, the accuracy obtained for each parameter, in particular 
the interesting one p c , involves a series of accumulated errors. Instead of 
taking care of p c , let's turn our attention to Aq(t) for a while. Figure 2 
shows a plot of this value function of r. 

One can see that A vanishes for the particular value r* = 0.984786(11). 
For the geometry we adopted, this is the equivalent of the above quoted Pin- 
son's number Ej, i.e. the universal value of our spanning probability at 
p c in the thermodynamic limit. Recently, Cardy JU] tried to calculate this 
value for an annulus, i.e. an L\ x L 2 square lattice where periodic boundary 
condition is adopted only along the direction of L\. The quoted universal 
value corresponds to the critical spanning probability between the two free- 
boundary lines, and depends only on the ratio r = L2/L1. Unfortunately, 
within his theoretical approach, Cardy was forced to leave out spanning con- 
figurations which also wrap along the direction of L 1; the periodic boundary. 
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Figure 2: Expansion parameter Aq, equation (1), as a function of r. The error 
bars obtained from the Chi-square fitting appear only within the smaller scale 
of the inset, near the special point r* = 0.984786(11) where Aq vanishes. 

Our geometry with the complete L x L torus divided in two halves by 
the lines i and i + L/2 corresponds to the r = 1/2 Cardy's geometry counted 
twice, in parallel. Thus, our r* can be related to Cardy's number C by 
1 - t* = (1 - C) 2 , i.e. C = 0.876657(45). Compared with Cardy's exact 
value C — C wrap = 0.7569977963 jTU] we found the contribution C wrap = 
0.119659(45) for the spanning configurations which also wrap around the 
periodic direction. 

The critical occupation probability p c can be obtained from our expansion 
(1) by fixing any value for r. The best choice is r = r*, for which the leading 
finite-size term in (1) becomes 

£-1-1/* instead of L~ x l v . In this case, we get 
p c = 0.59274621(33), where the error bar is estimated by the same above- 
quoted procedure of dividing the whole data set into S — 10 sub-sets. For 
some unknown particular reason, the wrapping probabilities adopted in P 
works better yet, the leading finite-size term in (1) being L -2-1 ^. In a 
numerical evidence for that behaviour is given, by plotting tl(p c )—t* against 
L (assuming some previously determined value for p c ) and verifying a power- 
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law dependence with an exponent very close to —2. Indeed, for our spanning 
probability, the same plot gives an exponent very close to —1. 

4 Conclusion 

We propose the finite-size expansion (1) for spanning probabilities in perco- 
lation, where Pl(t) is defined in figure 1. 

This approach can be used to calculate various critical quantities. We 
applied it to a particular geometry within the site percolation problem on a 
L x L square lattice, considered as a torus. Taking two parallel lines distant 
L/2 from each other, our spanning probability counts configurations for which 
these two lines are connected. The universal value r* = 0.984786(11) and 
the critical occupation p c = 0.59274621(33) are obtained. 

As a by-product, we propose the universal value C = 0.876657(45) for 
the critical spanning probability within a L x L/2 annulus, i.e. a square 
lattice with periodic (free) boundaries along the direction of L {L/2). This 
probability corresponds to all configurations for which the frontiers separated 
by L/2 are connected. Cardy [TU] determined the exact figure C — C wrap = 
0.7569977963, where C wrap corresponds to spanning configurations which also 
wrap along the direction of L, discounted in his approach. 

The whole computer time for obtaining these data was 20 thousand hours, 
on a dozen computers, typically powered by an Athlon 1GHz processor. 
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